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HIGHLIGHTS 


•  We  devise  a  control  strategy  to  maximize  battery  life  during  discharge. 

•  We  model  a  VRFB  system  to  predict  the  peak  power. 

•  The  model  considers  voltage  and  pump  losses  when  calculating  peak  power. 

•  The  model  accounts  for  mass  transfer  effect  and  non-uniform  vanadium  concentration. 

•  Electrode  porosity,  temperature,  and  vanadium  concentration  affect  peak  power. 
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The  vanadium  redox  flow  battery  (VRFB)  is  a  promising  grid-scale  energy  storage  technology,  but  future 
widespread  commercialization  requires  a  considerable  reduction  in  capital  costs.  Determining  the 
appropriate  battery  size  for  the  intended  power  range  can  help  minimize  the  amount  of  materials 
needed,  thereby  reducing  capital  costs.  A  physics-based  model  is  an  essential  tool  for  predicting  the 
power  range  of  large  scale  VRFB  systems  to  aid  in  the  design  optimization  process.  This  paper  presents  a 
modeling  framework  that  accounts  for  the  effects  of  flow  rate  on  the  pumping  losses,  local  mass  transfer 
rate,  and  nonuniform  vanadium  concentration  in  the  cell.  The  resulting  low-order  model  captures  bat¬ 
tery  performance  accurately  even  at  high  power  densities  and  remains  computationally  practical  for 
stack-level  optimization  and  control  purposes.  We  first  use  the  model  to  devise  an  optimal  control 
strategy  that  maximizes  battery  life  during  discharge.  Assuming  optimal  control  is  implemented,  we 
then  determine  the  upper  efficiency  limits  of  a  given  VRFB  system  and  compare  the  net  power  and 
associated  overpotential  and  pumping  losses  at  different  operating  points.  We  also  investigate  the  effects 
of  varying  the  electrode  porosity,  stack  temperature,  and  total  vanadium  concentration  on  the  peak 
power. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

With  the  rapid  increase  of  grid-connected  intermittent  renew¬ 
able  sources  such  as  wind  and  solar  energy,  it  is  becoming  crucial  to 
have  a  reliable  and  efficient  means  of  electrical  energy  storage  that 
can  maintain  stable  power  output  and  high  power  quality.  Among 
existing  technologies,  the  vanadium  redox  flow  battery  (VRFB)  is 
one  of  the  most  promising  solutions  for  grid-level  energy  storage 
because  of  its  fast  response  rate  and  high  efficiency  [1],  Further¬ 
more,  VRFBs  are  modular  so  they  can  easily  scale  to  meet  a  wide 
range  of  power  and  capacity  requirements  by  increasing  the 
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number  of  cells  in  the  stack  and  electrolyte  volume,  respectively. 
Still,  VRFBs  are  too  expensive  to  be  commercially  viable  [2], 

Operating  costs  can  be  reduced  by  improving  the  performance 
and  efficiency  of  VRFBs.  To  achieve  this  objective,  researchers  are 
designing  better  materials  [3]  and  cell  architectures  [4]  as  well  as 
developing  physics-based  VRFB  models  with  varying  levels  of 
complexity.  Shah  et  al.  [5—8  developed  first  principle  models, 
incorporating  fundamental  conservation  laws  and  electrode  ki¬ 
netics,  capable  of  predicting  2-D  distributions  of  concentration, 
overpotential,  and  current  density  in  a  unit  cell  for  a  range  of 
operating  conditions.  Xu  et  al.  [9]  developed  a  3-D  model  for 
investigating  different  flow  field  designs  to  reduce  overpotential 
losses  in  a  single  cell.  These  high  fidelity  models  are  valuable  for 
improving  our  understanding  of  the  important  physical  phenom¬ 
ena  that  affect  cell  performance  and  efficiency.  Nonetheless,  Rear 
et  al.  [10]  emphasized  the  additional  need  for  practical  stack-level 
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models  that  can  help  in  optimizing  the  design  of  large-scale 
systems. 

This  need  was  further  validated  through  an  economic  study  by 
Hittinger  et  al.  [11  ,  who  showed  that  improvements  in  capital  cost 
and  power  limit,  rather  than  efficiency  and  operating  cost,  have  the 
greatest  potential  impact  on  the  profitability  of  modular  battery 
storage  technologies  for  load  following  and  frequency  regulation 
applications.  Indeed,  optimal  battery  sizing  can  help  minimize  the 
amount  of  materials  needed  for  the  intended  power  requirement, 
resulting  in  lower  capital  costs.  For  this  reason,  there  is  a  strong 
incentive  to  pursue  research  not  only  in  engineering  new  materials 
to  operate  at  higher  power  densities,  but  also  in  developing 
modeling  and  simulation  tools  to  predict  the  peak  power  of  VRFB 
systems  so  that  their  full  operating  range  can  be  utilized  12],  Ma 
et  al.  [13]  conducted  experiments  to  characterize  the  performance 
of  a  specific  VRFB  system  and  devised  a  flow  rate  strategy  to 
improve  system  efficiency  at  various  states  of  charge  (SOC).  Xiong 
et  al.  [14]  took  a  more  general  approach  by  developing  a  lumped 
parameter  model,  accounting  for  temperature  effects  and  pumping 
losses,  to  predict  the  optimal  flow  rate  that  maximizes  system  ef¬ 
ficiency  under  different  operating  conditions  and  SOC.  Although 
this  model  is  computationally  practical  for  stack-level  optimization 
and  control  purposes,  it  lacks  sufficient  detail  in  the  voltage  model 
to  predict  battery  performance  accurately  when  operating  at  high 
power  densities. 

In  this  paper,  we  present  a  modeling  framework  that  accounts 
for  the  effects  of  flow  rate  on  the  pumping  losses,  local  mass 
transfer  rate,  and  nonuniform  vanadium  concentration  in  the  cell.  A 
sufficiently  high  flow  rate  is  desirable  for  increasing  the  local  mass 
transfer  rate,  thereby  reducing  the  concentration  overpotential. 
However,  an  optimal  flow  rate  strategy  should  consider  the  trade¬ 
off  between  reducing  overpotential  and  using  excessive  pump  po¬ 
wer.  A  quantitative  assessment  is  provided  in  this  paper  to  better 
understand  the  trade-off  between  these  two  competing  objectives 


in  order  to  come  up  with  an  optimal  flow  rate  strategy.  Depending 
on  the  flow  rate  and  reaction  rate,  a  lumped  parameter  model  may 
give  an  unsatisfactory  approximation  of  the  bulk  vanadium  con¬ 
centration  inside  the  cell.  To  improve  the  accuracy  of  the  lumped 
model  with  minimal  increase  in  computational  cost,  we  assume  a 
linear  profile  for  the  vanadium  concentration  along  cell  length  and 
calculate  an  average  open  circuit  voltage  and  overpotential.  The 
resulting  low-order  model  includes  sufficient  detail  to  capture  the 
internal  battery  losses  at  high  power  densities  while  remaining 
practical  for  stack-level  optimization  and  control  purposes.  First  we 
use  our  model  to  devise  an  optimal  control  strategy  that  maximizes 
battery  life  during  discharge  over  a  range  of  SOC  and  power  de¬ 
mand.  Assuming  optimal  control  is  implemented,  we  then  deter¬ 
mine  the  upper  efficiency  limits  of  the  given  VRFB  system  design. 
The  net  power  and  associated  overpotential  and  pumping  losses  are 
compared  quantitatively  at  different  operating  points.  To  further 
extend  the  analysis,  we  investigate  the  effects  of  varying  important 
design  parameters  including  the  electrode  porosity,  stack  temper¬ 
ature,  and  total  vanadium  concentration  on  the  peak  power.  During 
the  initial  design  stage,  these  insights  can  help  size  a  VRFB  system 
appropriately  for  the  intended  power  range  while  minimizing 
capital  costs. 


2.  Methodology 

We  consider  a  VRFB  system  as  shown  in  Fig.  1.  The  net  power 
output  of  a  VRFB  system  is  the  sum  of  the  stack  power  minus  the 
power  consumed  by  auxiliary  components  such  as  the  pump,  po¬ 
wer  electronics,  thermal  management  system,  and  control  system. 
In  our  analysis,  we  only  consider  the  internal  battery  overpotential 
and  pumping  losses  whereas  power  consumed  by  all  other  auxil¬ 
iary  components  is  neglected.  Pump  power  has  been  shown  in 
several  pilot  experiments  to  account  for  a  majority  of  the  parasitic 


Fig.  1.  Schematic  of  a  VRFB  system. 
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losses  in  a  VRFB  system  [1,9,13,16],  Thus,  the  system  net  power  can 
be  calculated  by 

Pnet  =  iE  oc  —  iE  loss  —  Ppump  (1) 


where  E0hm  is  the  ohmic  overpotential  and  tj pos  and  7jneg  are  the 
activation  overpotentials  at  the  positive  and  negative  half-cells. 

The  open  circuit  voltage  can  be  calculated  using  the  Nernst 
equation  as 


where  i  is  the  current,  Eoc  is  the  open  circuit  voltage,  £ioss  is  the 
internal  battery  overpotential  loss,  and  PpUmp  is  the  pump  power. 
The  following  section  explains  the  derivation  of  each  term  in  Eq.  (1 ). 


2.1,  Battery  model 


We  first  consider  the  lumped  parameter  model  of  a  unit  cell 
VRFB  system  developed  by  Shah  et  al.  [17  with  two  major  modi¬ 
fications.  First,  we  account  for  the  effect  of  flow  rate  on  the  local 
mass  transfer  rate  of  vanadium  ions,  which  could  cause  a  significant 
concentration  overpotential  at  high  current  densities.  Second,  we 
use  a  linear  profile  to  approximate  the  1-D  vanadium  concentration 
distribution  along  the  cell  length  and  calculate  an  average  open 
circuit  voltage  and  activation  overpotential  based  on  evenly  spaced 
grid  points. 

Dynamic  equations  can  be  obtained  by  partitioning  the  unit  cell 
VRFB  system  into  four  control  volumes  and  performing  a  mass 
balance  for  each  vanadium  species.  The  control  volumes  include 
the  positive  tank,  negative  tank,  positive  electrode,  and  negative 
electrode.  During  operation,  electrolyte  is  circulated  between  the 
storage  tank  and  the  electrode,  where  electrochemical  reactions 
occur.  The  electrochemical  reactions  that  occur  at  the  positive  and 
negative  electrodes  are  given  by 

Positive  electrode: 


t  ,  charge 

V02+  +  H20  ^  VO|  +  2H+  +  e- 

discharge 

Negative  electrode: 


(2) 


V 


■3+ 


charge 

discharge 


V' 


'2+ 


(3) 


Thus,  the  system  can  be  represented  by  eight  state  variables  in 
the  following  state  space  form: 

x(t)=Ax(t)Ui(t)  +  Bu2(t) 


A  = 


-  1  1  1  , 

eVP  4x4  eVP  4x4 


^tank 


'4x4 


^tank 


'4x4 


fi  =  4£L[-l  1  1  -1  0  0  0  0]r 

evpr 


x=\c2  c3  c4  c5  c2t  c3t  c4t  c5t  ] 

u  =  [q  nT 


where  E%os  and  £geg  are  the  standard  reduction  potentials  for  the 
reactions  at  the  positive  and  negative  half-cells,  R  is  the  molar  gas 
constant,  T  is  the  stack  temperature,  £  is  the  Faraday  constant,  and 
cH+  is  the  proton  concentration  in  the  positive  half-cell.  The  ohmic 
overpotential  includes  losses  in  the  membrane,  electrolyte  and 
current  collector  and  is  obtained  by 

^ohm  =  Em  +  Ee  +  Ec  (6) 

Each  ohmic  overpotential  term  is  calculated  respectively  by 


Em 


■  kVm 


(7) 


Ee  =j 


We 


(8) 


Ec  =j 


.  wc 


Oc 


(9) 


where  w  is  the  width  and  a  is  the  conductivity.  The  electrolyte 
conductivity  in  Eq.  (8)  includes  a  Bruggeman  correction  to  account 
for  the  electrode  porosity.  The  membrane  conductivity  varies  as  a 
function  of  temperature  according  to  the  Arrhenius  law  given  by 

[18] 


am  =  (0.5139A 


0.326)exp  1268 


1 

303 


(10) 


where  the  A  is  membrane  water  content.  The  membrane  is  assumed 
to  be  fully  saturated  with  A  =  22. 

The  activation  overpotential  including  mass  transfer  effects  can 
be  obtained  using  the  current-overpotential  equation  given  by  Ref. 

[19] 


(ID 


where  a  is  the  charge  transfer  coefficient.  The  exchange  current 
density,  jo,  is  calculated  by 
Positive  electrode: 

Jo, pos  =  kAmkpos  vT’4C5  (12) 

Negative  electrode: 


Here,  e  is  the  electrode  porosity,  Ve  is  the  electrode  volume,  Vtank 
is  the  tank  volume,  cvan  is  the  total  vanadium  concentration,  and 
/4x4  is  the  identity  matrix.  The  concentrations  of  V2+  and  V3+  in  the 
negative  half-cell  and  tank  are  represented  by  the  state  variables  c2, 
c2t,  C3,  and  C3t,  respectively.  The  concentrations  of  V02  and  V02+  in 
the  positive  half-cell  and  tank  are  represented  by  the  state  variables 
C4,  C40  c5,  and  cst,  respectively.  The  control  inputs  are  the  flow  rate, 
Q,  and  current  density,  j.  The  current  density  is  related  to  the  cur¬ 
rent  by  i  =  jAm,  where  Am  is  the  membrane  area. 

The  cell  voltage,  £ceii,  can  be  obtained  by 

£cell  =  Eoc  —  E0hin  —  V  pos  —  'in  eg  (4) 


Jo.neg  —  EAmknegx/C2C3  (13) 

The  reaction  rate  constants,  kpo s  and  kneg,  vary  with  temperature 
according  to  the  Arrhenius  law  given  by 


Epos  =  EpgS  exp  ( 


(E a 

1  1 

^ef“E 

(14) 


kneg  is  found  in  an  analogous  manner  to  Eq.  (14).  To  the  best  of 
the  authors'  knowledge,  only  the  positive  activation  energy  has 
been  reported  in  literature  [23]  so  it  is  assumed  that  the  negative 
activation  energy  takes  on  the  same  value  for  simulation  purposes. 
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The  cathodic  and  anodic  limiting  current  densities,  jiiC  and  j ka, 
are  the  maximum  current  density  limits  restricted  by  the  diffusion 
rate  of  vanadium  species  from  the  bulk  solution  to  the  electrode 
surface  and  are  related  to  the  mass  transfer  coefficient,  mvan,  by 
Positive  electrode: 


with  SOQank  at  the  cell  inlet  and  SOCcen  at  the  outlet  [22],  The  SOC 
at  a  distance  x  from  the  cell  inlet  can  be  obtained  by 

SOCcell(x)  =  SOCtank  -  (23) 

vicCvan  t-e 


jlc  — 

(15) 

Jl,a  =  — ^mmvanc4 

(16) 

Negative  electrode: 

Jl,c  =  ^mmvan^2 

(17) 

Jl,a  =  — ^mtnVanC3 

(18) 

All  vanadium  species  are  assumed  to  have  the  same  mass 
transfer  coefficient  [  15],  The  mass  transfer  coefficient  can  be  related 
empirically  to  the  electrolyte  velocity,  v,  by  Ref.  [20] 


mvan  =  2  x  10-V'4  (19) 

Assuming  that  the  positive  and  negative  sides  are  completely 
balanced  and  neglecting  vanadium  and  water  transfer  through  the 
membrane,  the  system  can  be  reduced  to  two  state  variables  by 
using  a  lumped  parameter  approximation  [17].  Accordingly,  the 
SOC  in  the  cell  and  tank  can  be  defined  respectively  as 


where  Le  is  the  electrode  length  in  the  flow  direction.  For  the 
simulations  in  this  paper,  we  use  a  resolution  of  N  =  5  equidistant 
grid  points  from  xi  =  0  to  X5  =  Le  to  calculate  an  average  open  circuit 
voltage  and  overpotential  by 


Eqc 


EPo°S 


~  Eoeg ' 


2RT 

~F~ 


+  N 


N 

Eln 

i=l 


f  SOCcen(Xi)  ^ 

Vl-SOCcell(x,Vj 


(24) 


and 


*  =  *!>  (25) 

i  l 

where  each  rp  is  a  function  of  SOCceii(x;)  and  satisfies  Eq.  (11). 

2.2.  Pump  power  model 

The  pump  power  required  to  operate  a  VRFB  system  is  calcu¬ 
lated  by 


SOCce|| 


C2  =  Q> 
Cyan  Cyan 


Cyan 


Cyan 


and 


SOCtank  =  C2t  =  C5t 

Cyan  LVan 


Cyan 


1 


c4t 

Cyan 


(20) 


(21) 


Eqs.  (5),  (11)— (13),  and  (15)— (18)  can  be  rewritten  in  terms  of 
the  two  new  state  variables,  SOCcen  and  SOCtank.  by  substitution  of 
Eqs.  (20)  and  (21).The  resulting  system  of  equations  then  becomes 


m  = 

=  Ax(t)Ui(t)  +  Bu2 

r  l 

i  i 

A  = 

eVe 

eVe 

1 

1 

.  e^tank 

£ ^tank  _ 

B  = 


Am 

eVeFcv3n 


1  T 


x  —  [  SOCcell  SOCtank  ] 

u  =  [Q_  j]T 


Because  the  tank  volume  is  typically  much  larger  than  the 
electrode  volume  in  a  utility-scale  VRFB  system,  SOCtank  changes 
much  more  slowly  than  SOCcen.  Assuming  that  SOCtank  remains 
constant  when  discharging  over  a  short  time  interval,  the  pseudo¬ 
steady  state  value  of  SOCcen  can  be  calculated  for  a  given  S0Ctank. 
current,  and  flow  rate  by  Ref.  [21] 


SOCcell  =  soctank  -  ^ - 

OiTCva  n 


(22) 


fpump  —  2??pumpQApt:ot 


(26) 


where  Aptot  is  the  total  pressure  loss,  7)pUmp  =  0.85  is  the  pump 
efficiency,  and  the  factor  of  two  accounts  for  the  two  pumps  needed 
to  circulate  the  positive  and  negative  electrolytes.  We  use  the  hy¬ 
draulic  circuit  configuration  based  on  a  typical  VRFB  system  design 
presented  by  Xiong  et  al.  [14]  to  determine  the  total  pressure  loss. 
In  this  section,  the  hydraulic  circuit  configuration  is  explained  and 
model  equations  are  summarized.  The  reader  should  refer  to  [14] 
for  more  specific  details  about  the  hydraulic  circuit. 

The  pressure  losses  in  a  VRFB  system  are  comprised  of  major 
and  minor  losses  that  occur  in  the  cell  stack  and  piping  system.  The 
major  loss  in  the  cell  stack  is  due  to  friction  in  the  porous  electrode 
and  is  calculated  by  using  a  modified  form  of  Darcy's  equation  [5,15] 


a  „  _  p,veLe 

^Fstack,  major  —  js 

where 


(27) 


I<  = 


dje3 


16k[<c(l  -  eY 


(28) 


Here,  p  is  the  electrolyte  dynamic  viscosity,  Le  is  the  electrode 
length  in  the  flow  direction,  Ae  =  eweLe  is  the  area  normal  to  the 
flow,  ve  =  QJAe  is  the  average  electrolyte  velocity  in  the  electrode,  d/ 
is  the  mean  fiber  diameter,  and  kKc  is  the  Kozeny-Carmen  constant. 
For  a  cell  stack  aligned  hydraulically  in  parallel,  the  pressure  drop 
across  each  electrode  is  the  same.  The  major  loss  in  the  piping 
system  is  due  to  friction  in  the  pipe  walls  and  is  calculated  by 


^Ppipe, major 


128pQf.p 


(29) 


Eq.  (22)  can  then  be  used  to  predict  the  overpotential  loss  as  a 
function  of  S0Ctank.  current,  and  flow  rate.  To  improve  the  accuracy 
of  the  lumped  parameter  model,  a  linear  profile  can  be  used  to 
approximate  the  1-D  SOC  distribution  along  the  length  of  the  cell 


where  Lp  is  the  pipe  length  and  Dp  is  the  pipe  diameter. 

Based  on  a  typical  flow  flame  design,  the  minor  losses  in  the  cell 
stack  can  be  approximated  by  summing  the  individual  minor  losses 
across  a  series  of  line  and  branch  flow  tees  14]: 
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(  v  1/  \ 

APstack, minor  =  2  X  nCells  (  ^lineP  +  ^branchP^  J  (30) 

Here,  J<iine  =  0.9  is  the  line  flow  tee  coefficient  and  /(branch  =  2  is  the 
branch  flow  tee  coefficient  for  threaded  fittings.  The  factor  of  two  is 
to  account  for  the  pressure  drop  at  the  stack  inlet  and  outlet.  nCeiis  is 
the  number  of  cells  in  the  stack.  In  a  typical  VRFB  piping  system 
[14],  the  minor  losses  can  be  approximated  by  summing  the  losses 
at  three  elbows  as  well  as  the  entrance  and  exit  of  the  tank: 

yp- 

^Ppipe, minor  =  3  X  /(C||K)W/'  ^  ^  x  ^"tnnk/;  2  (31) 

where  Keibow  =  1-5  is  the  elbow  friction  coefficient,  /(tank  =  0.97  is 
the  entry/exit  loss  coefficient  and  vp  =  Q/(7rD^/4)  is  the  average 
electrolyte  velocity  in  the  pipe. 

The  total  pressure  loss  in  the  VRFB  system  can  then  be  obtained 
by 


1  2  3  4  5  6  7 

Flow  Rate  (m3  s_1)  x  10 


(a) 


1400 


%  of  Max  Net  Power 

(b) 

Fig.  2.  (a)  Net  power  over  a  range  of  operating  points  (flow  rate  and  current)  when 
SOQank  is  30%.  (b)  Net  power  and  associated  pumping  and  overpotential  (OP)  losses  at 
three  different  operating  points:  Minimum  current  (left),  flow  rate  is  3.5  x  10  4  m3  s_1 
and  current  is  minimum  necessary  to  reach  given  net  power  (middle),  maximum 
current  (right). 


APtot  —  APstack,  major  +  ^Pstack, minor  ^Ppipe,  major  APpipe,  minor 

(32) 

and  substituted  into  Eq.  (26)  to  calculate  the  total  pump  power. 

3.  Results  and  discussion 

Fig.  2a  is  a  contour  plot  showing  the  net  power  as  a  function  of 
current  and  flow  rate  when  SOCtani(  is  30%.  The  pump  power 
calculation  is  based  on  a  20  cell  stack  and  the  net  power  output  is 
calculated  using  Eq.  (1).  The  parameter  values  used  for  the  simu¬ 
lations  are  provided  in  Table  1.  At  a  given  flow  rate,  the  limiting 
current  is  governed  by  the  rate  of  mass  transfer  of  vanadium  ions 
from  the  bulk  solution  to  the  electrode  surface.  As  the  current  ap¬ 
proaches  the  limiting  current,  the  concentration  overpotential  in¬ 
creases  significantly  causing  a  sharp  decrease  in  the  net  power. 
Consequently,  the  net  power  contour  is  cut  off  when  the  current 
exceeds  99%  of  the  limiting  current.  Fig.  2a  reveals  that  there  are  an 
infinite  number  of  possible  combinations  of  current  and  flow  rate  at 
a  given  SOCtank  that  yield  a  given  net  power  (except  the  maximum 
net  power,  for  which  only  one  operating  point  is  possible).  Physical 
intuition  suggests  that  the  optimal  operating  point  for  maximizing 
battery  life  during  discharge  uses  the  minimum  current  necessary 
to  provide  the  desired  net  power.  The  inherent  advantages  of 
operating  at  the  minimum  current  operating  point  (MCOP)  can  be 
realized  by  noticing  several  distinct  features  in  Fig.  2a.  Although  the 
MCOP  requires  slightly  more  pump  power  than  the  minimum  flow 
rate  operating  point  (MFOP)  at  a  given  net  power,  it  would  be  more 
practical  to  run  at  the  MCOP  because  the  MFOP  is  much  closer  to 
the  limiting  current.  Thus,  the  net  power  output  at  the  MFOP  would 
be  very  sensitive  to  changes  in  the  operating  point  due  to  con¬ 
centration  overpotential  effects.  The  plot  also  shows  that  the 
maximum  net  power  occurs  when  the  flow  rate  is  about 
3.5  x  1CT4  m3  s_1,  meaning  it  would  be  futile  to  increase  the  flow 
rate  further  beyond  this  point. 

Fig.  2b  compares  the  net  power  and  associated  battery  over¬ 
potential  and  pumping  losses  at  three  different  operating  points  on 
Fig.  2a  over  a  range  of  net  power  levels  when  SOCtani<  is  30%.  Each 
group  on  the  abscissa  represents  a  net  power  level  and  each  stack 
within  a  group  corresponds  to  a  different  operating  point  used  to 


Table  1 

Default  parameter  values  used  in  the  simulations. 


Parameter 

Description 

Value 

£o“ 

Positive  standard  reduction  potential 

1.004  V 

Eoeg 

Negative  standard  reduction  potential 

0.26  V 

R 

Gas  constant 

8.3145  J  mor1  K_1 

F 

Faraday  constant 

96845C  mor1 

Membrane  area 

0.01  m2  [17] 

Qf+ 

Proton  concentration 

4M  [17] 

wm 

Membrane  width 

1.25  x  10~4  m  |17] 

We 

Electrode  width 

4  x  10~3  m  [17] 

Wc 

Current  collector  width 

5  x  10~3  m  [17] 

<7e 

Electrolyte  conductivity 

100  Sm4  [17] 

<7C 

Current  collector  conductivity 

9.1  x  104  S  nr1  [17] 

a 

Charge  transfer  coefficient 

0.5 

href 

Positive  reaction  rate  @  Tref  =  293  K 

1.2  x  10-7ms-'  [23] 

href 

''•neg 

Negative  reaction  rate  @  Tref  =  298  K 

5.3  x  10~6  m  s  1  [24] 

Activation  energy 

52.4  x  103  J  mol-1  [23] 

Lp 

Pipe  length 

2  m  [14] 

Dp 

Pipe  diameter 

0.03  m  [14] 

P 

Electrolyte  density 

1400  kg  irr3  [14] 

P 

Electrolyte  dynamic  viscosity 

4.928  x  10~3  Pas  [15] 

Le 

Electrode  length 

0.1  m  [17] 

df 

Mean  fiber  diameter 

1.76  x  10~5  m  |15[ 

kra: 

Kozeny-Carmen  constant 

4.28  [15] 
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reach  that  net  power.  The  stack  on  the  left  corresponds  to  the 
MCOP.  The  stack  in  the  middle  corresponds  to  an  intermediate 
operating  point  where  the  flow  rate  is  3.5  x  1CT4  m3  s_1  and  the 
current  is  the  minimum  necessary  to  reach  the  given  net  power. 
The  stack  on  the  right  corresponds  to  the  maximum  current  oper¬ 
ating  point,  which  is  the  worst  possible  operating  point  during 
discharge.  This  plot  provides  an  illustrative  means  of  analyzing  the 
net  power  and  associated  losses  at  different  operating  points.  At  the 
MCOP,  the  flow  rate  is  just  high  enough  to  overcome  the  sharp 
increase  in  concentration  overpotential  without  excessive  pump 
power.  Consequently,  a  majority  of  the  losses  are  due  to  over¬ 
potential,  whereas  the  ratio  of  pumping  loss  to  net  power  is 
negligible  and  barely  visible  in  the  plot.  At  the  intermediate  oper¬ 
ating  point,  the  flow  rate  is  always  higher  than  necessary  to  reach 
the  desired  net  power,  resulting  in  a  greater  ratio  of  pumping  loss  to 
net  power.  The  maximum  current  operating  point  is  distinctly  the 
most  inefficient  of  the  three  operating  points  as  indicated  by  the 
significantly  higher  overpotential  and  pumping  losses. 


Flow  Rate  (m3  s  ’)  x  10 


(a) 


%  of  Max  Net  Power 


(b) 


Fig.  3.  (a)  Net  power  over  a  range  of  operating  points  (flow  rate  and  current)  when 
SOCtank  is  70%.  (b)  Net  power  and  associated  pumping  and  overpotential  (OP)  losses  at 
three  different  operating  points:  Minimum  current  (left),  flow  rate  is  3.5  x  10  4  m3  s  1 
and  current  is  minimum  necessary  to  reach  given  net  power  (middle),  maximum 
current  (right). 


Fig.  3a  and  b  are  for  the  case  when  SOCtank  is  70%  and  are 
analogous  to  Fig.  2a  and  b,  respectively.  When  SOCtank  is  70%,  Fig.  3b 
shows  that  the  ratio  of  pumping  loss  to  net  power  is  much  lower 
than  in  the  case  when  SOCtank  is  30%.  This  occurs  because  a  lower 
flow  rate  is  needed  to  reduce  the  concentration  overpotential  in  the 
cell  when  SOCtank  is  higher. 

In  the  discussion  so  far,  we  have  demonstrated  a  methodical  and 
physics-based  approach  to  devising  an  effective  control  strategy  for 
a  VRFB  system.  Assuming  that  an  optimal  control  strategy  was 
implemented,  it  was  then  straightforward  to  analyze  the  upper 
efficiency  limits  of  a  given  VRFB  system  design  by  comparing  the 
net  power  to  the  associated  overpotential  and  pumping  losses.  In 
the  next  part  of  this  section,  we  extend  the  analysis  by  showing 
how  the  maximum  net  power  of  a  VRFB  system  running  optimally 
is  affected  by  varying  important  design  parameters  including 
electrode  porosity,  stack  temperature,  and  vanadium  concentra¬ 
tion.  During  the  initial  design  stage,  this  insight  can  be  helpful  for 
sizing  the  VRFB  system  appropriately  to  handle  the  intended  power 
requirements  while  reducing  capital  costs. 

Fig.  4  shows  the  maximum  net  power  and  the  corresponding 
current  and  flow  rate  control  inputs  for  SOCtank  from  10  to  90%  with 
electrode  porosities  of  0.93,  0.8,  and  0.68.  This  range  of  electrode 
porosities  has  been  used  in  numerous  research  studies  and  pilot 
experiments  [5—8,15].  As  electrode  porosity  increases,  the 
maximum  net  power  increases.  This  trend  can  be  explained  by  Eq. 
(27),  which  predicts  that  head  loss  in  the  electrode  decreases  with 
increasing  porosity.  Therefore  with  higher  porosity,  the  pump 
consumes  less  power  at  a  given  flow  rate  and  the  limiting  current 
can  be  increased  further.  In  the  case  of  e  =  0.93  when  SOCtank  is 
high,  the  peak  power  is  achieved  at  a  slightly  lower  flow  rate 
compared  to  when  SOCtank  is  50%.  This  phenomenon  can  be 
explained  by  Eq.  (11),  which  predicts  that  the  minimum  activation 
overpotential  occurs  when  SOCceii  is  50%.  Using  a  lower  flow  rate 
not  only  decreases  the  pump  power,  but  also  makes  the  average 


Fig.  4.  Maximum  net  power  (top)  and  the  corresponding  current  (middle)  and  flow 
rate  (bottom)  control  inputs  for  SOCtank  from  10  to  90%  with  electrode  porosities  of 
0.93,  0.8,  and  0.68. 
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Fig.  5.  Maximum  net  power  (top)  and  the  corresponding  current  (middle)  and  flow 
rate  (bottom)  control  inputs  for  SOCank  from  10  to  90%  with  stack  temperatures  of 
313  K,  298  K,  and  283  K. 

SOC  in  the  cell  closer  to  50%  when  SOCtani;  is  high.  However,  the 
trade-off  is  that  decreasing  the  flow  rate  also  lowers  the  mass 
transfer  coefficient  and  therefore  the  limiting  current.  That  could 
likely  explain  why  the  same  trend  is  not  apparent  in  the  case  of 
e  =  0.8  or  0.68,  where  the  peak  power  is  achieved  at  a  relatively 
constant  flow  rate  when  SOCtank  is  high. 


Fig.  6.  Maximum  net  power  (top)  and  the  corresponding  current  (middle)  and  flow 
rate  (bottom)  control  inputs  for  SOCtank  from  10  to  90%  with  total  vanadium  concen¬ 
trations  of  2  M,  1.6  M,  and  1.2  M. 


Fig.  5  shows  the  maximum  net  power  and  the  corresponding 
current  and  flow  rate  control  inputs  for  SOCtank  from  10  to  90%  with 
stack  temperatures  of  318  K,  298  K,  and  278  K.  In  all  cases,  the 
electrode  porosity  is  0.93  and  the  total  vanadium  concentration  is 
2  M.  At  higher  temperatures,  the  increase  in  maximum  net  power 
can  be  attributed  to  the  decrease  in  activation  overpotential  ac¬ 
cording  to  Eqs.  (11)  and  (14).  Although  increasing  the  stack  tem¬ 
perature  increases  the  peak  power,  laboratory  studies  report  that 
vanadium  begins  to  precipitate  beyond  40  C  [14j.  These  results  can 
be  helpful  in  evaluating  the  necessity  and  determining  the  required 
capabilities  of  a  temperature  regulation  system  based  on  the 
climate  of  the  deployment  location. 

Fig.  6  shows  the  maximum  net  power  and  the  corresponding 
current  and  flow  rate  control  inputs  for  SOCtank  from  10  to  90%  with 
total  vanadium  concentrations  of  2  M,  1.6  M,  and  1.2  M.  In  all  cases, 
the  electrode  porosity  is  0.93  and  the  stack  temperature  is  298  K. 
The  maximum  net  power  evidently  increases  with  a  higher  total 
vanadium  concentration.  However,  the  flow  rate  required  to  ach¬ 
ieve  maximum  power  decreases  when  SOCtank  is  above  50%  for  2  M 
vanadium  concentration,  whereas  the  flow  rate  remains  relatively 
constant  for  the  1.6  M  and  1.2  M  cases.  With  a  lower  total  vanadium 
concentration,  a  higher  flow  rate  is  required  to  keep  the  average  cell 
overpotential  low.  Depending  on  the  total  vanadium  concentration 
being  used,  it  can  be  important  to  consider  the  flow  rate  strategy 
when  trying  to  take  full  advantage  of  the  peak  power  when  SOCtank 
is  high. 

4.  Conclusions 

In  this  paper,  we  developed  a  physics-based  model  incorpo¬ 
rating  battery  overpotential  and  pumping  losses  that  can  be  used  as 
a  valuable  tool  for  optimizing  the  design  and  control  of  a  stack  level 
VRFB  system.  Sufficient  detail  was  included  in  the  overpotential 
model  to  account  for  the  effects  of  flow  rate  on  the  limiting  current 
and  on  the  non-uniform  vanadium  concentration  in  the  cell. 
Therefore,  system  performance  can  be  predicted  accurately  even  at 
high  current  densities.  First  we  used  the  model  to  devise  a  control 
strategy  that  maximizes  battery  life  during  discharge.  Assuming  the 
optimal  control  strategy  was  implemented,  we  then  determined 
the  upper  efficiency  limit  of  a  given  VRFB  system  at  different  SOC 
and  net  power  levels.  The  net  power  and  associated  overpotential 
and  pumping  losses  were  also  compared  quantitatively  at  different 
operating  points.  Furthermore,  we  investigated  the  effects  of 
varying  the  electrode  porosity,  stack  temperature,  and  total  vana¬ 
dium  concentration  on  the  peak  power  output  for  SOCtank  from  10 
to  90%.  For  a  certain  combination  of  VRFB  system  design  parame¬ 
ters,  we  found  that  a  non-intuitive  flow  rate  strategy  should  be 
followed  to  take  full  advantage  of  the  peak  power  available  when 
SOCtank  is  high.  The  methodology  and  analysis  demonstrated  in  this 
paper  can  readily  be  adapted  to  optimize  other  critical  parameters 
in  the  design  process  and  help  to  further  reduce  capital  costs. 
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